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ABSTRACT 

The quasiparticle excitation spectrum of an isolated vortex in a clean layered d-wave superconductor is cal- 
culated, using both the Bogoliubov-de Gennes and the quasiclassical Eilenberger equations. Analytical solutions 
are obtained within the model of step- variation of the gap function, adjusting the Caroli, de Gennes and Matricon 
approach for low- lying excitations in cuprates. A large peak in the density of states in the "pancake" vortex core 
is found, caused by the two-dimensional and strong coupling nature of high-temperature superconductivity. 
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1 INTRODUCTION 

Enhancing interest has been focused on quasiparticle spectra properties around vortices in both conventional 
and unconventional superconductors. Density of localized states near the Fermi surface in the vortex core in 
conventional type-II superconduntaip is equivalent to that of the normal metal area with radius of the order of 
the coherence length £ ~ 10 nrnJilcm High-temperature superconductors (HTS) are layered, having a cylindrical 
Fermi surface, the superconductivity is of d-wave and strong coupling type, the coherence length being much 
smaller, £ ~ lnm. This implies a crucial difference in the quasiparticle spectra in the vortex core between 
conventional superconductors and HTS. 

Numerous computations of the quasiparticle spectra around two-dimensional vortices have been performed 
recently. For s-wave pairing, Rainer et al. used the quasiclassical theorya and Hayashi et al. self-consistentljz. 
solved the Bogoliubov-de Gennes equations in the quantum limit.cl For d-w&ve pairing, Schopohl and MakiQ 
and Ichioka et alE analyzed the vortex structure within the quasiclassical theory, and Morita et alE solved 
the Bogoliubov-de Gennes equations, neglecting the "quantum effects", i.e. non-commutativity between the 
quasiparticle momentum and coordinate, obtaining the four-fold symmetry of the vortex core. A fully self- 
consistent -numerical approach, including " quantum effects" , has been done by Morita et alu and Franz and 
TesanovicEJ The common feature of layered superconductors is the appearance of a large peak in the quasiparticle 
density of states (DOS) in the vortex core. rT,his is in agreement with the scanning tuapeling microscopy (STM) 
observations on NbSe2, done by Hess et <^~f^ and on YBCO by Maggio-Aprile et oL.EJ However, in BSCCO the 
gap structure was found by Renner et aZ. J13 without any signature of the quasiparticle bound states in the vortex 
core. Possible explanation is either the existence of the pseudogap ip .the normal metal phase,£3 or the low-lying 
states are not localized but extended along the gap node directions. OllJ 

Since the Andreev bound states can transport charge currents, unlike the bound states in a potential well, 
supercurrents can flow through the vortex, strongly influencing its dynamics a .T.n .particular, this could be very 
important for understanding anomalous electric and thermal Hall conductivity!^'! 16 !! 17 ! and the quantum magnetic 
dc Haas- van Alphen oscillations £3 

In this paper, both the Bogoliubov-de Gennes and the Eilenberger quasiclassical equations are solved 
analytically, within the model of step- variation of the gap function, connecting the bound states in the vortex 
core and in the normal metal cylinder embedded in the same superconductor. Since the bound states spectrum is 
not quasicontinuous, for 2D and strong coupling superconductivity, this approach modifies the Caroli, de Gennes 
and Matricon results for low-lying states. Our results indicate that DOS has either one large maximum or a gap 
structure, depending on the value of fcj?£. 



2 BOGOLIUBOV DE GENNES EQUATIONS 

An efficient method for calculating local spectral properties is based on the Bogoliubov-de Gennes equations 

{ — (-ihV - -A] - E F \ u + Av = Eu, 
{ 2m V c / J 

-|^- (-ihV + -E F ^v + A*u = Ev, (1) 

where u and v are the Bogoliubov amplitudes and A is the gap function. 

For a normal-metal cylinder of radius r c , embedded in a layered superconductor along the z-axis and in 
the absence of magnetic field, A = 0, in cylindrical coordinates (r, (p, z), z \\ c, we assume the gap function in the 
form 

0, r < r c 



A = A M ) = j A( , )5 r> ; c , (2) 
where 

A(9) = A cos(26») (3) 

for d-wave pairing. The self-consistency condition and non-cpmmutativity between coordinate and momentum 
components are neglected, 9 being the polar angle in k-space.B In this case, Eqs. ([!]) can be rewritten in the form 
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with n = 0,±1,±2, .... 

For r < r c , Eq. (0) has the following real- valued solution, finite at r = 



(6) 



f ± = A ± J n (q ± r) . (7) 



Here, J„ is the Bessel function of the first kind, A± are arbitrary real numbers, and q± = y kp ± 2mE/U 2 . 
For r > r c , the solution can be chosen in the form 



f±=J2 c ±h> (8) 



x 

where A and f\ are the eigenvalues and the corresponding eigenfunctions of L, 

Lfx = A/a. (9) 
For bound states with energies below the superconducting gap, \E\ < |A(0)| 



A = ±i^A 2 (9)-E 2 , (10) 
and f\ are the solutions of the Bessel equation 
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with ql = {2m/h 2 )y/A 2 {6) - E 2 . Introducing 

k = ^k 2 F + iq 2 (fcfc > 0), 
the real- valued solutions, which tend to zero as r — > oo, are 
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Continuity of /+ , /_ and their first derivatives at r — r c yields the system of algebraic equations for A± , 
SRB and 3_B, which has a nontrivial solution if 
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Here £ = ^A 2 (0) -£ 2 , J„(±) = J n (q±r c ), HP = H ( n\kr c ), J'{x) = dJ(x)/dx, H'{x) = dH(x)/dx. Bound 
state energies are calculated numerically from Eq. (Eq). 



For an isolated "pancake" vortex in HTS in a low magnetic field, H ~ H c i <C H C 2, applied along the z-axis, 
we assume the same model for A(r, 6), Eq. (|2|), taking AkO and the gauge 



A = A(r, 0)e" 



(16) 



due to the flux quantization. The core radius r c is now r c (6), of the order of the BCS coherence length £o(#) = 
hv F /Tr\A(0)\. Taking 
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where 2v is an odd integer, Eq. ([|) approximately describes the vortex core, if L — > Z/ and — > where 
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describe the influence of the screening currents flow around the vortex on the low-lying bound states. In this way 
the solution for the vortex core is related to the solution for normal metal cylinders with (9-dependent r c . This 
approximation is close to the approach in Ref. [1], adjusted for HTS, where only one or two lowest bound states 
appear. 

Density of bound states 

N b (E) = ^ f^deJ^HE-Eiie)), (20) 
2W Y 

where i is n for a cylinder and v for a vortex, is calculated approximately for E > 0, taking for a cylinder 

£ n (0) w£„(O)|cos(20)|, (21) 



and for a vortex 



Therefore, for a cylinder 
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and for a vortex 
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where a = E' u (0), and f3 — \is\h 2 /2mt;(0) 2 . Low-lying bound states are well described by taking the constant 
energy shift _Ln Eq. (|22|), whereas for larger v the shift becomes zero, E'(6) w E(9), as shown by perturbative 
calculations o 



3 EILENBERGER EQUATIONS 



In this section, we consid^jij the same model in the quasiclassical approach. We start with the Eilenberger 
equations in the following forn 



2e 

2huj n + hw F • ( V — 
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(25) 
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Here g = gn(r,vp,uj n ) and / = f^(r,vp,oj n ) represent the normal and the anomalous Green function respec- 
tively, A = A(r,vp) is the gap function, oj n — irksTfin + 1), n = 0, ±1,±2..., are the Matsubara frequencies, 
and vp is the Fermi velocity vector. The function ft is defined by ft, (r, vp,u> n ) — /*-|-(r, — v>,u; n ). A and / are 
connected by the self-consistency equation. 

For a homogeneous and isotropic superconductor, solutions of the Eilenberger equations are 



A f A* hu> n 



(27) 



where 



e 2 = |A| 2 + (hu n ) 2 . (28) 

For a pancake vortex in (r,tp) plane, using the model gap function, Eqs. (||) and (|l6|), and denoting the 
coordinate along vp by s, and along £_x vp by p (Fig.l.), in the gauge with real-valued gap function, Eqs. ( |25| ) 
and (p6|) can be rewritten in the formtil 
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(30) 
(31) 



where r 2 = p 2 + s 2 . 



Fig.l. Trajectory passing at distance p from the vortex center. 

Normal metal cylinder in zero magnetic field is described by Eqs. (p9|)-(^f|) without ip/r 2 term. For r < r c , 

f = Fe- K " s , g = G, (32) 
where kq = 2u> n /vp. For r > r c and k — 2e n /hvp 

/=</)+ <f ie - KS , g = (g)+T 1 e- KS ) for s > 0, (33) 
/ = (/)+ $ 2 e KS , g=(g)+T 2 e KS , for s < 0. (34) 

Eqs. ©-© imply 
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From the continuity condition for / and g at ±sq = ± \Ac — > f° r r < r c the normal Green function is 



G = 



?Lw„ cosh(ftoso) + e « sinh(KoSo) 
S.w n sinh(KoSo) + £«, cosh(KgSo) 



(36) 



For a vortex, approximating p/r 2 by p/r 2 (9) 1 the solution of Eq. (p9|)-(|3l|) can be obtained from Eq. (f3(i|), 
by changing w n -» w B ', 

cosh(K s ) + e' n sinh(rtoSo) ^ 



G 



where 



hu)' n sinh(KQSo) + e' n cosh(K so) ' 



(38) 



n " 2r c 2 (0)' 
and k' = 2u>' n /vF- 

Performing an analytical continuation of G by Tiuj n —* — zi? + 0, E being the quasiparticle energy with 
respect to the Fermi level, the retarded propagator g R (E,p, 9) is obtained. 
For a normal-metal cylinder, the bound states spectrum is given by 
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For a vortex and low- lying states, El Eq. p» holds for E — > E 1 , 
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In terms of reduced variables E/A Q E, y / A 2 (9) - E 2 /A -> s, ^ E 2 - A 2 (9)/A -> e, p/r c -> p, 
2so/<o -> s , the angle-resolved partial DOS (PDOS) is obtained from N(E,p,9) = %tg R (E,p,9). 



For the normal-metal cylinder 



N(E,p,0)/N(O) = 9(e 2 ) ^ , + 

e 2 cos 2 (Es )+ E 2 sm (Es ) 

+ (e 2 ) 7r| ^ (g)l S (Esm(Es ) - e cos(Es )) , (41) 

where (5 is the Dirac function, 9 is the step-function, and N(0) — m/Zirfi 2 is the normal metal D9S for one spin 
orientation. Averaging PDOS over the cylindrical Fermi surface gives 
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Finally, after spatial averaging over the cylinder area 7rr 2 , DOS is 



MS) = - / N(E,p)y / l-p 2 dp. 
t Jo 



(43) 



For the vortex, in Eq. (|l|) .E £ + E sign£, £ = ^|p 1/ |w jF /2A o r 2 (0), p v = v/k F , and Eqs @ and @ 
become more complicated. For small-radius vortices in HTS, kp^o(0) ~ 1, where only one trajectory through the 
vortex core is allowed, with p — pi = 1/2/cf, the quasiclassical DOS, with one large peak, is qualitatively the 
same as for the normal-metal cylinder. 



4 RESULTS AND DISCUSSION 

Our results are illustrated in Figs. 2-6. for A /E F = 0.4, corresponding to k F £ (0) = 2E F /nA = 1.5, 
which is a reasonable value for cuprates, providing one distinguishable peak in DOS in the vortex core. 

For a normal-metal cylinder of radius r c — 5£o(0) embedded in a layered d-wave superconductor, seven 
distinguishable localized states are found (Fig. 2). This number obviously decreases with approaching the gap 
function nodes in the k-space. For a cylinder of radius r c = £o(0), the minimal vortex core radius, only two bound 
states are distinguishable, as shown in Fig. 3. Note that the lowest level of the quasicontinuum bound states 
calculated from quasiclassical theory (dashed curve in Fig. 3) is higher than in the quantum approach. However, 
for a d-wave pancake vortex, only one localized state can be seen (Fig. 4.) in the core of radius r c (9) = £o(0). 
Note that only this state (y = 1/2) is shifted, according to the arguments from the more rigorous perturbation 
theory approach. E3 

Density of bound states, calculated within the quantum model, for the normal-metal cylinder and for the 
vortex core, is shown in Fig. 5. Only one peak in DOS appears (except for the usual discontinuity at the gap 
edge) for the chosen value k F ^ (0) = 1-5. In the quasiclassical approach the total DOS, including the unbound 
states, is calculated for the normal-metal cylinder for the same parameters and shown in Fig. 6. Qualitatively, 
this is the same as the dashed curve in Fig. 5, despite the fact that the quantization of bound states is not taken 
into account. Note that in the quasiclassical approach the peak energy is systematically shifted towards higher 
energies. 

In conclusion, the approximate analytical method is developed in the same manner as the Caroli, de Gennes 
and Matricon approach, adjusted for the case of small-radius pancake vortices, where the quasiparticle spectrum 
is not quasicontinuous. Two-dimensional and strong-coupling nature of superconductivity in cuprates implies 
that the density of localized quasiparticle states inj-tlie. vortex core is rather different from the normal metal DOS. 
Comparing our results with the experimental dataj!3ll3 it seems that only one large peak in DOS is distinguishable 
for 0.7 < k F £ (0) < 1.5, as in YBCO, and none (gap structure) for fcp£o(0) < 0.7, as in BSCCO, in consistency 
with a larger gap for the latter .E3 This can be the simplest explanation of the experimental observations. 



Fig. 2. Quasiparticle energy E/Aq vs the momentum orientation angle in a normal-metal cylinder, calculated 
from Eq. @, for fc F £o(0) = 1.5 and r c = 5£ (0). Shaded region: unbound states. 



Fig 3. Quasiparticle energy E/Aq vs the momentum orientation angle in a normal metal cylinder, calculated 
from Eq. ([l5|), for kp£,o(0) = 1.5 and r c = £o(0). Dashed curve: the lowest bound state (p = 0), calculated from 
Eq. ( p9[) within the quasiclassical theory. Shaded region: unbound states. 



Fig. 4. Quasiparticle energy E/Aq vs the momentum orientation angle in the vortex core of radius r c — ($(6), 
calculated from Eq. jl5|), for fc.F£o(0) = 1.5. Shaded region: unbound states. 



Fig. 5. Density of bound states in the vortex core r c = £,o(9) (solid curve) and in the normal-metal cylinder 
r c = £o(0) (dashed curve), calculated from Eqs. ( |23| ) and (p4|), for /c,f£o(0) = 1-5. 



Fig. 6. Quasiparticle DOS in the normal metal cylinder, calculated from Eq. (^3|), for r c — £o(0) and &f£o(0) 

I ..*). within the quasiclassical theory 
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